Phase-slip avalanches in the superflow of ^He through arrays of nanopores 
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Recent experiments by Sato et al. jj have explored the dynamics of '*He superflow through an ar- 
ray of nanopores. These experiments have found that, as the temperature is lowered, phase-slippage 
in the pores changes its character, from synchronous to asynchronous. Inspired by these experi- 
ments, we construct a model to address the characteristics of phase-slippage in superflow through 
nanopore arrays. We focus on the low-temperature regime, in which the current-phase relation for 
a single pore is linear, and thermal fluctuations may be neglected. Our model incorporates two 
basic ingredients: (1) each pore has its own random value of critical velocity (due, e.g., to atomic- 
scale imperfections), and (2) an effective inter-pore coupling, mediated through the bulk superfluid. 
The inter-pore coupling tends to cause neighbours of a pore that has already phase-slipped also to 
phase-slip; this process may cascade, creating an avalanche of synchronously slipping phases. As the 
temperature is lowered, the distribution of critical velocities is expected to effectively broaden, owing 
to the reduction in the superfluid healing length, leading to a loss of synchronicity in phase-slippage. 
Furthermore, we find that competition between the strength of the disorder in the critical velocities 
and the strength of the inter-pore interaction leads to a phase transition between non-avalanching 
and avalanching regimes of phase-slippage. 

PACS numbers: 67.40.-w, 67.40.Hf, 83.60.Df 



Quantum phenomena at the macroscopic scale have been 
exhibited in a variety of physical settings, including su- 
perconductors 0, the hehum superfluids and, more 
recently, Bose-Einstein condensates (BEC) in atomic 
vapours 0|, and perhaps even supersolids Q. Spectac- 
ular examples include persistent currents in multiply- 
connected superconductors and superfluids, and the 
Josephson effects JpL first realised in superconducting 
tunnel junctions [3, l3] and, later, between two we akly 
coupled BECsJi, as weh as in superfluid ^He [l3[ll[l2 
and ^He [H [l|| weak links. 

Widely-used descriptions of these systems are based 
upon the unifying concept of spontaneous symmetry 
breaking implemented within the framework of the 
Landau theory of phase transitions. In this framework, 
one introduces a suitable macroscopic Ginzburg-Landau 
order parameter to characterise the broken symmetry 
state. One major virtue of order-parameter based de- 
scriptions is the view they afford of important topologi- 
cal processes, such as vortex nucleation and growth, and 
the mechanisms for dissipation that ensue. 

Superconducting macroscopic quantum phenomena 
(MQP) have been employed in many technologies, in 
fields ranging from biotechnology to radio astronomy. 
Primarily based on the DC Josephson effect — the es- 
sential element in superconducting quantum interfer- 
ence devices (SQUIDs) — these include SQUID amplifiers, 
metrology, biology and medicine, nondestructive testing, 
geomagnetism, magnetic measurements and microscopy, 
and prospecting for oil and minerals 16]. 

The superfluid analogue of the SQUID differs in an es- 
sential way from its superconducting counterpart: it is 
sensitive to interference amongst (charge-neutral) atoms 
rather than (charged) electrons. As a result, it has con- 
siderable potential for use in measuring angular veloci- 



ties rather than magnetic fields a nd, therefore, for use as 
an ultra-sensitive gyroscope [rallj. As with SQUIDs, 
the essential element is the DC Josephson effect, which 
has been realised in "^He 0, ^| and, more recently, in 
4He [a Ei. The latter is expected to be more use- 
ful, technologically, owing to its relatively high operating 
temperature, which exceeds that of ^He-based systems 
by a factor of two thousand. 

Recent experiments by the Berkeley group have ex- 
plored a system that comprises two superfluid reservoirs 
coupled via an array of pores, over a wide range of tem- 
peratures [23. It was observed that, for progressively 
lower temperatures below the critical temperature T\ for 
the transition to the superfluid state, the system passes 
through several regimes. Just below Tx, for sufficiently 
small supercurrents, the array was observed to function 
as a single effective Josephson junction having an essen- 
tially sinusoidal current-phase relation. For lower tem- 
peratures, below roughly a few millikelvin below T\, it 
was observed that the single-pore dynamics becomes ir- 
reversible. This is due to the dissipative phase-slips that 
occur whenever the superflow velocity in a pore reaches 
a critical value [23. 123L 123. |25| , which we believe to be 
specific to that pore. At the high-temperature end of 
this regime, there is a narrow interval of temperatures, 
of width roughly 10 mK, within which all pores appear 
to phase-slip simultaneously, and which Sato et al. 
refer to as the synchronous regime. At lower tempera- 
tures, down to roughly 160 mK below T\, it appears that 
simultaneity in phase-slipping is lost; Sato et al. refer to 
this as the asynchronous regime. It is believed that tech- 
nologically viable devices can be designed to function in 
both the Josephson (i.e. reversible) and phase-slippage 
(i.e. irreversible) regimes, provided a large enough frac- 
tion of pores slip sufficiently simultaneously to produce 
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FIG. 1: Schematic of the model system. Left: the location 
of the pore array on the membrane is indicated by the black 
region, and the phases on hemispheres at infinity are labelled 
and . Centre: slice through the membrane, with pores 
being represented by breaks in the membrane (white). Right: 
boundary conditions on hemispheres at the openings of the 
i-th pore. 



a measurable "whistle" at the Josephson frequency |2(il |. 
We remark that, in the setting of multi-link supercon- 
ducting devices, it has been shown that the irreversible 
regime can be utilised for magnetic-field and related 
phase-sensitive measurements |27l |. 

In this Paper, we address the issue of the transition 
from synchronous to asynchronous phase-slip dynamics 
of an array of nanopores connecting a pair of supcrfluid 
reservoirs. The main ingredients of our description are 
pores that have random, temperature-dependent criti- 
cal velocities, along with an effective inter-pore coupling 
mediated via superflow in the reservoirs. We develop a 
model that incorporates these ingredients, and analyse it 
via both a mean-field approximation and exact numerical 
analysis for arrays consisting of a relatively small number 
of pores. Thus, we identify two effects: (a) strong dis- 
order washes out the synchronicity of phase slips, which 
leads to the loss of the "whistle" ; and (b) if the disorder 
is sufficiently weak, the phase-slip dynamics undergoes 
a disorder-driven phase transition between avalanching 
and non-avalanching regimes (see Fig. ^ . This enables 
us to obtain the current-phase relation, as well as the 
amplitude of the current oscillations at the Josephson 
frequency corresponding to a fixed chemical-potential dif- 
ference. We believe that this model and our analysis of it 
captures the essential physics taking place in the Berke- 
ley group's experiments 1]. 

Disparate physical systems featuring competition be- 
tween quenched disorder and interactions, such as sliding 
tectonic plates 28] , random-field magnets , and solids 
with disorder-pinned charge-density waves (CDWs) js^, 
l3l| . are well known to show phenomena analogous to 
the avalanching-to-non-avalanching transition described 
here. The model and analysis that we employ are similar 
to those that have been applied in the aforementioned 
settings [s^l- 

Basic model — The system we wish to describe consists 
of two reservoirs of supcrfluid ^He, separated by a rigid 
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FIG. 2: Phase diagram showing avalanching and non- 
avalanching regimes of the phase-slip dynamics, as a function 
of the effective pore strength J/27rro and disorder strength 
w. The diagram was computed via our mean field theory, 
in the large-array limit, with top-hat critical phase twist (or, 
equivalently, the velocity) distribution of width w. Inset: The 
fraction of pores that have phase slipped, as a function of the 
control-parameter ($^ — il>^)/27r. Comparison between nu- 
merics for a 25 X 25 array of pores (points) and our mean-field 
theory (lines) with Gaussian distribution of critical velocities. 

barrier, embedded in which is an array of pores, as shown 
in Fig. n] We shall specialise to the case of an x iV 
array of pores, each having radius ro, situated at the sites 
of a square lattice of lattice parameter £. It is straight- 
forward to extend our analysis to other array geometries. 
It is convenient to regard the system as comprising three 
components: two are bulk components [i.e., the left (L) 
and right (R) reservoirs]; the third consists of the su- 
pcrfluid inside the pores. We describe the state of the 
bulk helium in terms of the supcrfluid order-parameter 
phase fields x'^/^ (r). In doing this we are neglecting ef- 
fects of amplitude excitations of the order parameter, in- 
cluding vortices. In contrast, within the pores we retain 
both amplitude and phase degrees of freedom. We imag- 
ine controlling the system by specifying the two phases, 
^^/^^ to which the phases in the left and right reser- 
voirs tend, far from the pore array. We believe that 
this level of description allows us to capture the following 
important elements: (a) pores that exhibit narrow-wire- 
like meta-stable states, these states being connected by 
phases slips; and (b) interactions mediated through the 
bulk superfiuid in the two reservoirs, which couple pairs 
of pores to one another and also couple the pores to the 
control phases, <I>'"/^. 

To summarise, we describe the bulk supcrfluid helium 
reservoirs by the phase-only Hamiltonians 

HUn^Ksf rf3^|v;,WR.(r)|^ (i) 

^ Jh/R 

where Ks = h^ns/m is the supcrfluid stiffness, in which 
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rig is the superfluid number density and m is the mass 
of a ■*He atom. To account for phase-shppage processes 
within a pore, which arise from vortex hues crossing the 
pore, we shall use a modified phase-only model that keeps 
track of the number of phase slips. Therefore, we take 
the energy of the superfluid inside the i*^ pore to be 



H, 



-Jl 



27rn,- 



(2) 



in which J = nr'^/d accounts for the geometry of the 
pore, where d is of the order of the membrane thickness, 
and is the phase of the order parameter in the vicin- 
ity of the left/right side of the i"^ pore. Furthermore, rii 
counts the net number of phase slips that would occur in 
the i**^ pore if the system were to progress from a refer- 
ence state in which the phase were uniform throughout. 
This description of the pores must be supplemented by 
specifying the critical velocities at which phase slips oc- 
cur, as we shall discuss in detail below. 

To reduce the model to one involving only the phase- 
differences across the pores and the phases at infinity, we 
minimise the energy Q in the reservoirs, which forces 
T^L/R^j.^ to obey the Laplace equation. Focusing on the 
left reservoir, the appropriate boundary conditions are 
(see Fig.^: (a) the phase on the hemisphere at infinity 
is <i>'^; (b) no current flows through the membrane surface 
between the pores, i.e. Vj.x'^ = there; and (c) the phase 
on the hemisphere of radius rg centred at the opening of 
the i^^ pore is specifled to be . The choice of surface 
for the last of the boundary conditions is made for con- 
venience, as it simplifles the resulting mixed boundary 
value problem whilst enforcing the physical condition of 
continuity of the phase in the vicinity of the pore open- 
ing. 

To solve this mixed boundary value problem we appeal 
to its electrostatics analogy, in which the phase and the 
superfluid stiffness Kg respectively play the roles of the 
scalar potential and the permittivity eg. To compute the 
energy at the minimum in terms of the boundary data 
we apply the divergence theorem to Eq. to obtain the 
energy in terms of a surface integral: 



K, 

2 -idh 



f dS-x^Vx^-:^ / dVx^VV, (3) 



where the volume (i.e. last) term vanishes because at the 
minimum-energy conflguration obeys V^x'" = 0. To 
simplify the evaluation of the surface integral we make 
use of the fact that the energy is invariant under global 
shifts of the potential. Thus, we lower all potentials by 
, which eliminates the contribution from the (left) sur- 
face at inflnity. What remains are the contributions from 
the membrane: those from between the pores give zero, 
owing to the zero-flux boundary condition; those from 
those from the hemispherical surfaces covering the pores 
give 



(4) 



where, by assumption, the phase x^ takes the constant 
values (f)\ on the hemispherical surfaces, and the remain- 
ing surface integrals have been replaced (via the analogue 
of Gauss's law) by half of the (analogue of) the charge 
enclosed. For pore i this charge is qi = Kg J dS ■ Vx'", 
where the integral is taken over the spherical surface that 
completes the hemispherical one. Reflecting the prob- 
lem across the membrane's left surface, we see that we 
are looking for half of the electrostatic energy of a set of 
suitably charged metal spheres at potentials — and 
what remains is to determine the charges Qi. To do this, 
we consider the i^^ sphere: in the limit ^ 3> tq (i.e. ignor- 
ing di- and higher-order charge- multipoles), the potential 
on its surface obeys 
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(5a) 
(5b) 



where \rij\ is the distance between the i"^ and j^^ pore 
(sphere) centres, and Cy is the analogue of the capaci- 
tance matrix. By solving Eq. (|5a|) for the g^'s and elim- 
inating them from Eq. Q), we arrive at the combined 
energy of the left and right superfluid reservoirs and the 
superfluid in the pores: 



where we have, without loss of generality, restricted our 
attention to spatially anti-symmetric states, for which 

We complete the description of our model by spec- 
ifying the single-pore dynamics, and thus the mecha- 
nism by which energy is dissipated in the pores. The 
superfluid velocity Vi in a pore of thickness d is de- 
fined by the phases at the pore openings: Vi = 
hV4>i/m Ri h{2(t>\ — /dm. Correspondingly the 

current through the pore is given by 



li = -J— (20, 



27rn,- 



(7) 



When the velocity through the i pore exceeds its critical 
value Vc^i (or, equivalcntly, <j>f — nni exceeds <j>c.i), ^ vor- 
tex line nucleates and moves across the pore, which de- 
creases the phase-difference across the pore by 2tt. Thus, 
immediately after a phase slip, the velocity through the 
^th pQj.g jg decreased by Av = 2TTh/md, a quantum of 
superfluid velocity drop for a pore having flxed phases 
at the openings, i.e. Vi Vi — Av. However, after a 
very short time, controlled by the speed of sound, the 
system balances the superflow in the bulk reservoirs and 
through the various pores. Therefore, after this relax- 
ation process is complete, the actual drop in the super- 
fluid velocity through the i"^ pore is always less than Av. 
To flnd the configuration of the superfluid after a phase 
slip, we note that the phase-difference along a path from 
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the left hemispherical surface at infinity through the i 
pore to the right hemispherical surface at infinity drops 
by 2tt, whilst the phase-difference along a path through 
any other pore remains unaffected. In our model, we im- 
plement this kind of phase-slip event by sending Ui to 
Tii + l (assuming all flow is to the left) and finding a new 
set of values for all of the (pf^s by minimising the total 
energy, Eq. ©. 

Not too near T\ the experimentally-observed 
temperature-dependence of the critical velocity for su- 
perflow through a pore is approximately linear l34| : 
Vc{T) = Wc(0)(l - T/To), where Tq ss 2.5 K. This 
T-dependence is most closely reproduced by the ther- 
mal nucleation of half-ring vortices mechanism |35| . 
However, the linear T-dependence does not hold in the 
temperature regime that we are primarily interested 
in; instead, the critical velocity is proportional to the 
superfluid stiffness there jSHSSEll: 

Va{T)^h/m^{T)^K,{T), (8) 

where ^(r) ~ ^o(l - T/Tx)^^/^ is the temperature- 
dependent healing length. 

Extrinsic effects are known to reduce critical velocities 
from the intrinsic values discussed above. We hypoth- 
esise that in the Berkeley group's experiments the ex- 
trinsic effects originate in atomic-scale roughness of the 
pore walls, and play a pivotal role in generating critical- 
velocity variability amongst the pores. This variability 
is expected to be temperature dependent because only 
roughness on length-scales longer than ^(T) can substan- 
tially perturb the order parameter and thus influence the 
critical velocities. Hence, at higher temperatures the im- 
pact of surface roughness is expected to be weaker and, 
correspondingly, the distribution of critical velocities is 
expected to be narrower. Thus, lowering the temperature 
has the important effective consequence of increasing the 
effective disorder [svj . 

Implications of the model — We shall work at fixed (neg- 
ative) difference Afj, in the chemical potential between 
the reservoirs, so that the control parameter evolves 
linearly in time, according to the Josephson-Anderson 
relation 

$L^„$R^„^i. (9) 

As $L_(j)R gj-Q-^vs^ so do the superfluid velocities through 
the various pores, punctuated at regular intervals by ve- 
locity drops associated with the phase-slip processes. As 
the total energy of the state is periodic in <i>^ with period 
TT, the total current through the array must be a periodic 
function of time with the period given by the Josephson 
frequency wj — l^^/h. Due to the randomness of the 
critical velocities amongst the pores, the velocities in the 
various pores do not reach their critical values simulta- 
neously. Therefore, the pores having the smaller critical 
velocities (i.e. the weaker pores) tend to slip first. If the 
distribution of critical velocities is sufhciently narrow, the 
array may, as we demonstrate below, suffer an avalanche. 



1 




'8 10 12 14 16 

cOj t / 2re 



FIG. 3: Traces of total current through an array of pores as a 
function of time, at various disorder strengths l38l . The traces 
were computed at fixed chemical potential difference A/^, via 
our mean-field theory, with Gaussian distributions of critical 
phase-twists (jia of widths g and mean = Stt. The corre- 
sponding critical current, in the absence of disorder, is given 
by Ic ~ 2KsJ4>c/h. As the disorder strength is increased, the 
amplitude of current oscillations decreases. The sharp drops 
in the current, which correspond to avalanches, disappear for 
cr/<^c > 0.052. 



By an avalanche we mean that when the weaker pores 
slip, superflow through the neighbouring pores that have 
yet to slip increases, due to the inter-pore interaction, 
and this drives them to their own Vc,ij causing a cas- 
cade of phase slips in which an appreciable fraction of 
pores in the array slip. Experimentally, avalanches are 
reflected in a periodic series of sharp drops in the total 
current through the array of pores as a function of time. 
Time-traces of the total current in the avalanching and 
non-avalanching regimes are contrasted in Fig. |21 

For arrays having a small number of pores, the quasi- 
static state of (mechanical) equilibrium may be numer- 
ically determined, as the control parameter — 
evolves parametrically. As a consequence of the long- 
range nature of the inter-pore couplings, the array dy- 
namics is well approximated by a mean-field theory. We 
shall describe a numerical approach first, and then the 
mean-field theory. 

The numerics take as input: J, the Wc.i's, and the (non- 
inverted) matrix C~j^ [see Eq. (j5b|l ] . which itself depends 

on ro, Z, and N. At each time-step, is incremented, 
and the new 0["'s are obtained. If the superflow through 
any pore is found to now exceed its critical velocity, that 
pore phase-slips (i.e. its value of rn is incremented by plus 
unity); next, the various ^^'s are recomputed using the 
new set of n^'s, and the program goes back to recheck 
if any other pore now exceeds its critical velocity. This 
continues until no new phase-slips are found to occur, 
at which point the control parameter is incremented and 
the procedure is repeated. 

Next, we construct a mean- field theory. We assert that 
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FIG. 4: Left- and right-hand sides of the averaged form of 
Eq. nut . The circles indicate the two main solutions: as the 
control-phase is increased, the RHS curve moves up and to 
the left, such that there is a discontinuous jump from one 
solution to the other corresponding to an avalanche. 
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where Q{(j)c) is the probability distribution of half-critical 
phase twists in the pores 0c = dmvc/fi- The averaged 
form of Eq. 1101 can be solved graphically, by plotting 
the left- and right-hand sides as functions of {(f>^); see 
Fig. 01 It is evident from this graphical approach that 
whenever the maximum slope of the right-hand side, 
2TrJd{n)/d{(f>^) fails to exceed the slope of the left-hand 
side, {A + 2J)((/)^), the self-consistency condition yields 
a unique solution for the average phase {(j)^), and that 
this phase evolves continuously with the (increasing) con- 
trol phase This corresponds to the non-avalanching 
regime. By contrast, whenever the maximum slope of 
the right-hand side does exceed the slope of the left- 
hand side, the self-consistency condition no longer yields 
a unique solution for {(j)^). Instead, as the control-phase 
increases, the continuous evolution of {(f)^) is punctuated 
by jumps, which occur when pairs of solutions merge and 
disappear. These jumps reflect avalanching behaviour, 
and we refer to this as the avalanching regime. 

One can use this mean-field theory to construct a 
phase diagram that demarcates avalanching and non- 
avalanching regimes, for any choice of disorder distribu- 
tion Q. For the case of a top- hat distribution of critical 
phase-twists, a simple inequality defines the avalanching 
regime: 

2ttJB , , 

" - "'^ = (A-f2J)(C + 2J)' 



the control phase-difference is monotonically increasing 
in time, so that phase differences c^y^ — (j)^ are always in- 
creasing and the superflow in the pores undergoes only 
n.i-increasing phase slips. Then we may proceed by se- 
lecting an arbitrary pore z, and minimising the energy © 
with respect to <j>}' ; by replacing the '/'^(^i) 's by the mean- 
field value {(j)^) we obtain an equation for the phase at 
the i^^ pore: 

4>]:{C + 2J) - B{(f)^) = A^^ + 2TTJni, (10) 

where A = J2i Qj, B ee Cu - A, and C = Cu 39]. (The 
necessary inversion of C^^ can readily be accomplished 
either analytically, by transforming to Fourier space, or 
numerically.) By averaging the left- and right-hand sides 
of this equation over sites, we arrive at a self-consistency 
condition on (c/)^). Next, by assuming self- averaging with 
respect to the disorder in the critical velocities, we may 
replace the average over sites by an average over disor- 
der. This procedure amounts to replacing, in the above 
condition, (/)^ by {(p^) and by 



(11) 
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where w is the width of the top hat, and Wc is its critical 
value; such a phase boundary is exemplified in Fig. |21 
In experiments, one can explore this phase diagram by 
tuning the temperature, which, as we have hypothesised 
above, effectively tunes the strength of the disorder. 

Thermal fluctuations of the phases wash out the 
disorder-driven phase transition if they exceed the width 
of the disorder distribution. To avoid this, the temper- 
ature has to be smaller than the energy cost of winding 
the phase of a single pore by the critical disorder width, 
i.e. < Ks{T) Cwl/2. We estimate that this inequal- 
ity is satisfied provided that T is not too close to the A- 
point, i.e. Ta - T > 3mK in the setting of Ref. [3. Fur- 
thermore, we have assumed that the current-phase rela- 
tion is linear, which is true provided that T\ — T> 5 mK, 
as was measured for a similar setup in Ref. [l^ . 

To test the results of the mean-field theory, we have 
compared them to results from a numerical investigation 
performed on a finite, periodic lattice [39j |. For various 
widths of the disorder distribution (which we have taken 
to be Gaussian), the fractions of pores that have phase- 
slipped, as a function of control-parameter, are shown 
in the inset of Fig. |21 As predicted by the mean-field 
theory, avalanches occur when the distribution of criti- 
cal velocities is narrow but not when it is broad. At a 
critical strength of the disorder, which separates these 
two regimes, the discontinuity in the mean-field fraction 
of slipped pores just vanishes. Moreover, at this critical 



e 



A^^ + B{(l)^) -nkC 
C + 2J 



e 



A'l>^ + B{(l)^)~Tr{k-l)C 
C + 2J 



6 



0.05 r 




ity, Eq. (|SJ, to include disorder, by assuming that 

h 



VcAT) 



0.02 0.04 0.06 0.08 0.1 
Relative disorder, a/(|) 



0.12 



FIG. 5: Amplitude of the oscillations of the current A/, and 
drop in the current caused by an avalanche, as functions of 
the disorder strength. A and N-A indicate the avalanching 
and non-avalanching regimes, respectively. Inset: Solid line: 
amplitude of the current oscillations as a function of temper- 
ature, using the disorder model described in the text. Dashed 
hne: expected amplitude in the absence of disorder. 



disorder there are always values of the control phase at 
which the response of the system diverges. The mean- 
field and numerical results appear to agree with one an- 
other rather well, as one can see from Fig. |21 at least in 
the vicinity of the critical point. 

The two main results of our Paper are summarised in 
Fig. [SI The blue curve shows the dependence of the am- 
plitude of the current oscillation A/ (i.e. half the distance 
between smallest and largest current during a single pe- 
riod in Fig. 131) on the disorder strength. As the disorder 
becomes stronger, the phase-slips in the various pores 
become more asynchronous, and the oscillations in the 
current gradually disappear. The red curve shows the 
dependence of the current drop caused by an avalanche 
(i.e. the height of the vertical drop in current in Fig. 
on the disorder strength. The current-drop plays the 
role of an order-parameter in a second-order phase tran- 
sition tuned by the disorder strength. As the disorder be- 
comes stronger, the order-parameter decreases, becoming 
zero at the critical disorder strength corresponding to the 
transition from the avalanching to the non-avalanching 
regime (i.e. ctc ~ 0.0520c)- 

Comparison with experiments — In their experiments 0, 
the Berkeley group have measured the amplitude of the 
"whistle" (i.e. the amplitude of the current oscillations 
A/) as a function of temperature at a fixed chemical 
potential difference. The experiments found that at Tx 
there are no oscillations of the current. As the temper- 
ature is lowered below Tx, AI first increases, and then 
gradually decreases. To obtain the behaviour of AI as a 
function of T we extend the model for the critical veloc- 



m(^(T)-f xO' 



(13) 



where Xi is the temperature-independent characteristic 
scale of the surface roughness in the i*'' pore, which we 
take to have a Gaussian distribution For Tx — T > 
5 mK, we can compare the results of our modelling to 
those of the experiments. The general features are repro- 
duced: the initial increase in AI is associated with an 
increase in the superfluid fraction; the gradual decrease 
at lower temperatures is due to the loss of synchronicity 
amongst the pores, which is caused by the effective in- 
crease in the strength of disorder. To demonstrate these 
features, we plot the amplitude of current oscillations as 
a function of temperature; see the inset in Fig. [S] for 
which we have set Sx = A nm. 

We also note that the general features of the current- 
ws.-time traces. Fig. El are similar to those of the type 
III experiments described in Ref. [ij. In both, as the 
temperature is lowered (i.e. the disorder is increased), 
the avalanche gradually disappears, and then so do the 
oscillations in the current. 

Concluding remarks — Motivated by recent experiments 
performed by the Berkeley group on superflow through 
nanopore arrays 0, we have developed a model to de- 
scribe phase-slip dynamics of such systems. The main 
features of our model are effective inter-pore couplings, 
mediated through the bulk superfluid, as well as random- 
ness in the critical velocities of the pores, the latter being 
effectively controlled through the temperature. 

Within our model, we find that the competition be- 
tween (a) site-disorder in the critical velocities and (b) ef- 
fective inter-pore coupling leads to the emergence of 
rich collective dynamics, including a transition between 
avalanching and non-avalanching regimes of the phase- 
slip dynamics. We identify a line of critical disorder- 
strengths in the phase diagram, at which there is a di- 
vergent susceptibility, in the sense that near to this line 
small changes in the control parameter can lead to large 
changes in the fraction of phase-slipped pores. 

Our model reproduces the key physical features of 
the Berkeley group's experiments including a high- 
temperature synchronous regime, a low-temperature 
asynchronous regime, and a transition between the two. 
We therefore feel that the model captures the essential 
physics explored in these experiments. 

In a forthcoming paper we shall extend our approach in 
order to address the transition from the Josephson regime 
to the avalanching phase-slippage regime described here, 
by including thermal fluctuations together with a gener- 
alised description of pore energies that is valid in both 
regimes. 
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